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It is well known that magnification bias, the modulation of galaxy or quasar source counts 
by gravitational lensing, can change the observed angular correlation function. We investigate 
magnification-induced changes to the shape of the observed correlation function w(8), and the an- 
gular power spectrum Ce, paying special attention to the matter-radiation equality peak and the 
baryon wiggles. Lensing effectively mixes the correlation function of the source galaxies with that 
of the matter correlation at the lower redshifts of the lenses distorting the observed correlation 
function. We quantify how the lensing corrections depend on the width of the selection function, 
the galaxy bias b, and the number count slope s. The lensing correction increases with redshift and 
larger corrections are present for sources with steep number count slopes and/or broad redshift dis- 
tributions. The most drastic changes to Ci occur for measurements at high redshifts (z > 1.5) and 
low multipole moment (I < 100). For the source distributions we consider, magnification bias can 
shift the location of the matter-radiation equality scale by 1-6% at z ~ 1.5 and by z ~ 3.5 the shift 
can be as large as 30%. The baryon bump in 6 2 w(6) is shifted by < 1% and the width is typically 
increased by ~ 10%. Shifts of > 0.5% and broadening > 20% occur only for very broad selection 
functions and/or galaxies with (5s — 2)/b > 2. However, near the baryon bump the magnification 
correction is not constant but is a gently varying function which depends on the source population. 
Depending on how the w(9) data is fitted, this correction may need to be accounted for when using 
the baryon acoustic scale for precision cosmology. 

PACS numbers: 98.80.k,95.30.Sf 



I. INTRODUCTION 

The galaxy two-point angular correlation function and 
its spherical harmonic transform, the angular power spec- 
trum, provide information about dark matter clustering. 
These statistics have scale dependent features which may 
be used as "standard rulers" to measure cosmological dis- 
tances 0, 0. Features in the three-dimensional power 
spectrum which appear at a comoving wave number k 
will appear in the angular power spectrum at redshift zq 
at multipole £ ~ fcx(zo), where x( z o) is the comoving 
distance to z$. One such feature is the peak in the power 
spectrum that separates the modes which entered the 
horizon during radiation domination from those modes 
that entered during matter domination. This scale is de- 
termined by the horizon size at matter-radiation equal- 
ity. Also present in the power spectrum are the so-called 
baryon acoustic oscillations (BAO). This series of wig- 
gles in Fourier space is a signature of acoustic oscilla- 
tions in the photon-baryon fluid that was present in the 
early universe. The location of the peaks of the BAO 
depends on the horizon size at the time of recombination 
[IB H, IE 0i H> • This scale is robustly measured from 
the cosmic microwave background [ToL UA\ , and thus can 
serve as a "standard ruler"; given the scale of the BAO 
a measurement of the angular size of the BAO at some 
redshift will yield the angular diameter distance to that 



redshift . In fact the first measurements of this peak have 
occurred recently [H, [13, [TH, [H, In recent years, 
much effort has gone towards using these features in the 
correlation function for precision measurements fl2l |. 

Gravitational lensing changes the observed number 
density of galaxy or quasar sources - an effect called mag- 
nification bias [ill, [l^, H3| ■ (Hereafter, the terms 'galaxy' 
and 'quasar' can be considered synonymous.) It is well 
known that magnification bias modifies the gala xy angu- 
lar correlation function [U [H, HI H3, HEIM iM 111- 
In this paper we extend the previous analyses by inves- 
tigating and quantifying how magnification bias changes 
the shape of the angular two point correlation function, 
and its spherical harmonic transform the angular power 
spectrum, paying special attention to important features 
such as the turnover in the power spectrum and the BAO. 

Corrections from gravitational lensing enter as follows. 
The two-point function is measured from the galaxy num- 
ber density fluctuation. 



The effect of gravitational lensing is to alter the area of 
the patch of sky being observed and to change the ob- 
served flux of the source. Both effects change the mea- 
sured galaxy number density, the first by changing the 
area, the second by changing the number of sources ob- 
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served in a flux limited survey. Together these effects are 
called magnification bias [Hf, [2(| Hl[ . To first order 
(i.e. the weak lensing limit), they lead to a correction 
term 5^ being added to the intrinsic galaxy fluctuation 



Sn = S„ + Sn 



(2) 



With this term the observed autocorrelation function be- 
comes, 



(SnSn) = (S g 5 g ) + (SgS^) + (5M + (S^Sn 



(3) 



The galaxy-galaxy term (S g S g } depends on the matter 
distribution at the source galaxies. The magnification 
terms, especially (5^5^), depend on the matter distri- 
bution spanning the range between the sources and the 
observer. 

The lensing of high redshift quasars by low redshift 
galaxies have been detected confirming the presence of 
magnification bias for these systems. The most recent 
measurements of this effect are discussed in [25|, |H, |29| 
(see also [13] )• Discussions of earlier measurements can 
be found in the references therein. These measurements 
work by cross-correlating the angular positions of galax- 
ies/quasars at widely separated redshifts. Here, we will 
focus on the angular correlation of gala xies at similar 
redshifts. As was first pointed out by [30(, magnification 
bias alters observations of the 3D clustering of galaxies. 
In two separate papers [3l|, [HJ , we further consider the 
effect of magnification bias on the 3D correlation function 
and power spectrum, which turn out to have qualitatively 
new and interesting features. 

While this paper was in preparation, a paper ad- 
dressing weak gravitational lensing corrections to baryon 
acoustic oscillations was posted [331 ] . In that paper they 
consider the effect of magnification bias (a first order cor- 
rection) and stochastic deflection (a second order correc- 
tion) on BAO in the real space correlation function for a 
delta function source distribution. In contrast, we con- 
sider the effect of magnification bias on the angular cor- 
relation function and the angular power spectrum for an 
extended source distribution. As we will show the width 
of the source distribution strongly affects the magnitude 
of the magnification bias correction. While some con- 
clusions we reach are addressed in [33[, the observables 
considered and the analyses here are different. 

This paper is organized as follows. In section |TT] we 
present expressions for the angular auto-correlation func- 
tion and the angular power spectrum when magnification 
bias is included. In section HTT1 we discuss the factors af- 
fecting the relative magnitude of the magnification bias 
terms compared with the galaxy term. In section IIVI wc 
examine the effect of magnification on the shape of the 
angular power spectrum for a few different source distri- 
butions. In section [V] we consider the effect of magni- 
fication on the baryon bump in the angular correlation 
function. In I VII we conclude and discuss the implications 
of our work for future surveys. 



II. ANISOTROPIES AND CORRELATION 
FUNCTIONS 

We consider the two-dimensional galaxy fluctuation in 
direction 9 integrated over a selection function with mean 
redshift zq 



5 n (9,z ) = 



n(9, z ) - n(z ) 



n{z ) 



(4) 



where n is the mean number of galaxies in the redshift 
bin. We include this label, zq, to remind the reader that 
the angular fluctuation depends on the source selection 
function. 

When the effect of magnification bias is considered, 
the expression for the net measured galaxy overdensity 
becomes a sum of two terms (Equation [2]). In the expres- 
sions that follow and throughout this paper we assume a 
flat universe, generalizing to an open or closed universe is 
fairly straightforward. The first term in Equation[2]is the 
intrinsic galaxy fluctuation integrated over a normalized 
selection function given by W(z, zq), 



S g {9,z a ) 



dzb(z)W(z,z )8(x(z)9,z). (5) 



Where x( z ) 1S the comoving distance to redshift z, S = 
G° — P)/P is th e matter overdensity and b(z) is a bias 
factor relating the galaxy number density fluctuation to 
the matter density fluctuation S g ~ b(z)S. For simplic- 
ity we have assumed that the bias is scale independent, 
this should be accurate for k < 0.05/iMpc _1 , however 
at smaller scales it may be important (see for example 
[35l]). If we further assume that b(z) is slowly varying 
across W(z, zq) then we can set b(z) — b(zo) and pull it 
out of the integral in Equation 

The second term in Equation [5] is the correction due 
to magnification bias, 

f°° c 

5 M (d,zo) = J dz-^yVi0( X (z)0,z o ) (6) 
x x (z) r dz' (5s(z') - 2) X{Z ' ] ~ A X(Z) W(z', zo). 



H{z) is the Hubble parameter and c is the speed of light 
[19| . The magnification bias term depends on the Lapla- 
cian of the gravitational potential 4> (with respect to the 
comoving coordinates in the direction perpendicular to 
9) and the slope of the number count function. For a 
survey with limiting magnitude m this is 



d\og w N{< m) 
dm 



(7) 



If s(z) is slowly varying across W(z, zq), then we can set 
s(z) = s(zq) and the expression for 5^ becomes 



6»{fl,zo) = (5*0%) -2) 



x / dz 



j(z,zo)V 2 ± <p(x(z)9,z ). (8) 
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Where we have introduced the lensing weight function 



dz' 



X(z') 



-W(z',z ). (9) 



The lensing weight function can be thought of as roughly 
proportional to the probability for sources in W(z, zq) 
to be lensed by density perturbations at z. The lensing 
weight function increases in magnitude as Zq increases 
and is peaked at a redshift z corresponding to about half 
the comoving distance to zq. 

On scales k ^> aH Poisson's equation can be used to 
relate the gravitational potential to the matter fluctua- 
tion 

-fc 2 0(k,z) = ^n m (i + z)S(k, z) (io) 

where k is the comoving wave vector j k is its magnitude 
and f2 m is the matter density today [3J|. . 

We consider the two-point correlation function of the 
galaxy overdensity 

w nn (6,z Q ) = (5 g (6,zo)6 g (O,zo))+2(5 g (0,z a )5^(0,zo)) 
+ (S^e\z )S f ,(e\z Q )) 

= w gg (9,z ) + 2w gil (9,z ) +w l _ ifl (9,z ) (11) 

where cos 9 = 9-6. In what follows we will consider the 
Legendre coefficients of the correlation functions. These 
are defined as 



Wgg(9,Z ) = ^ 



21+1 
47T 



C 99 (z )P e (cos9) (12) 



where Pf(cos#) are the Legendre polynomials. The Leg- 
endre components of the observed correlation function, 
Wnn(0, zq) will of course be a sum of all the terms, 

Cr(zo) = Cf(z ) + 2Cf (zo) + C^(z ) (13) 

We calculate the angular power spectra, C 99 , C 9 ^ and 
, using the Limber approximation 36] which is accu- 
rate for t > 10. To simplify the expressions and provide 
some insight into the effects of the magnification terms 
we introduce the following functions 



<P"(z,z ) = b(z a ) 2 W(z,z a ) 2 

r^z,z ) = 

9 
4 



(14) 
(15) 



(5s(zo) - 2) 2 ^ n _iL_ 5 2 (Zj Zf))(1 + z) 2 



9 ^(Z,Z O ) = y/(j)33{z,Zo)^(z,Zo) 



(16) 



With this notation the angular correlation functions 
can all be expressed as 



dZ H ^^z o) p("z) (17) 



o X(z) 2 c 



X(z)' 
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FIG. 1: The redshift distributions normalized by the galaxy 
bias (Eq. [14115)) . 0" /b 2 (dashed line), 20 9fl /6 2 (dot-dashed 
line) and M/J /b 2 (dotted line) shown here for Gaussian se- 
lection functions centered at zq = 1.5 and zo = 2.5 with 
a = 0.15. These distributions determine the angular power 
spectra through Equation 1171 The total auto-correlation 
C" n /b 2 will depend on (<j> sg + 2^ + c/>^)/b 2 (solid line). 



where xx symbolizes gg, gfi or fifj, and P{k = £/x{z), z) 
is the matter power spectrum. 

For illustration we use Gaussian selection functions 
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of various widths to demonstrate how the shape of the 
selection function alters the effect of magnification. 

We assume a flat ACDM cosmology with fl m = 0.27, 
= 0.73 and fl b = 0.0224//i 2 as the fractional en- 
ergy densities in matter, vacuum and baryons today. 
The Hubble constant Hq = lOO/i km/s/Mpc is set to 
h = 0.7, the fluctuation amplitude to as = 0.8, and 
scalar spectral index n s = 0.95. For the linear power 
spectrum we use the Eisenstein and Hu transfer func- 
tion [8j. In a few places (e.g. Figured]) we discuss the 
no-BAO spectrum, this is calculated using the BBKS 
transfer function [37| with a modified shape parameter 
T = ft m hexp[-Q b (l + V2h/Q m )} [38[. The non-linear 
evolution of both power spectra is calculated using the 
prescription of Smith et. al. (3{|. The nonlinear power 
spectrum is used in all plots and discussions. 

Unless otherwise stated we set t he g alaxy sample de- 
pendent ratio (5s — 2)/b = 1 (see [U Eoj on how this 
varies with galaxy sample and redshift). The first cor- 
rection term, C 9fl /b 2 is linear in this quantity while the 
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second term, Cf/b 2 is quadratic. Thus the magnitude 
of the net correction, calculated here cannot in gen- 
eral be scaled by (5s — 2)/b. However, at low redshifts 
2Cf/b 2 dominates over Cf/b 2 and at higher redshifts 



Cr > 2C. 



The redshift of the transition between 



the two regimes increases with the width of the selection 
function. Specifically when a = 0.07, 2Cf < Cf for 
z > 1.5, when a = 0.15, 2Cf < Cf for z > 2.0 and 
when a = 0.30, 2Cf < Cf for z > 2.5. 



III. AMPLITUDE OF THE LENSING 
CORRECTIONS 

The relative magnitude of the lensing magnification 
terms, Cf and Cf, to the intrinsic galaxy term, Cf , 
depends on several things: first the galaxy-sample depen- 
dent quantities b and s, and second the selection function 
and cosmological quantities in Equations [14] UH Here we 
discuss how these quantities affect the relative magni- 
tudes of Cf, Cf and Cf. 

The magnification terms are scaled by the galaxy- 
sample dependent factors (5s — 2)/b and (5s — 2) 2 /b 2 . 
Thus for a given galaxy bias, galaxies residing on the 
steep end of the luminosity function will have larger 
magnification corrections. However, if s < 2/5, the 
galaxy-magnification term Cf is negative, while the 
magnification-magnification term, Cf is always positive. 
Unless otherwise stated we fix (5s — 2)/b = 1 so both 
terms are positive. 

The difference between the galaxy-galaxy, galaxy- 
magnification and magnification-magnification terms can 
be seen more clearly by considering equations [T4UT71 
From Equation [171 we see that all the angular power 
spectra are simply integrals of the power spectrum over 
some redshift distribution <f) xx (z,zo). Examples of these 
distributions for a selection function with a = 0.15 are 
shown in Figure [TJ In the case of c/) 99 (z,zq) this dis- 
tribution is determined by the selection function. This 
suggests that one can think of the magnification terms as 
adding to Cf 3 measurements of the correlation function 
with "selection functions" determined by <j) 9li {z,zo) and 
^fz,z a ). 

Figure Q] shows that 4> 9fi (z,zo) is peaked not too far 
from where the selection function is peaked, and is sim- 
ilar in shape to the selection function. Thus we expect 
Cf to be similar in shape to Cf 9 . On the other hand, 
the magnification-magnification term 0^ is peaked at a 
much lower redshift than the selection function and is 
also much more broadly distributed in redshift. Thus 
we expect Cf to be peaked at lower £ (larger angular 
scales) than Cf 9 because it is probing structure at lower 
redshifts which is nearer to the observer and occupies a 
larger angle in the sky. Additionally, since cj)^^ (z , zo) is 
rather broadly distributed, we expect sharp features such 
as baryon oscillations to be smoothed out in Cf . 

One might infer from looking at Figure[T]that the mag- 
nification bias terms Cf and Cf will be quite small 




100 1000 
I 



FIG. 2: The angular power spectra divided by the no-BAO 
power spectra for a selection function of width a = 0.15. The 
dashed line is Cf alone (no magnification case), the solid line 

is Cl n = Cf + 2C| M + Cf. 



compared to Cf . However, the other quantities in the 
integrand of Equation[T7]decrease with redshift so the rel- 
ative magnitudes of <fi 99 , (j) 9 ^ and MM do not completely 
determine the relative magnitudes of Cf,Cf and Cf. 



IV. THE SHAPE OF THE ANGULAR POWER 
SPECTRUM 

Using equations [T4HT51 and [17] we calculate the an- 
gular power spectra using Gaussian selection functions 
centered at a variety of redshifts. Because the galaxy- 
magnification cross term Cf depends strongly on the 
width of the selection function, wc consider several 
different widths given by cr = 0.07,0.15,0.30. The 
magnification-magnification term Cf is largely inde- 
pendent of the width of the selection function, but the 
galaxy-magnification term Cf increases with increasing 
width of the selection function. The effect of increas- 
ing a is therefore to increase the overall contribution of 
magnification to Cf. 

In Figure [2] we show the angular power spectrum with 
and without magnification bias for a selection function 
of width a — 0.15. In this figure the power spectra are 
normalized by the no-BAO power spectrum, which is cal- 
culated neglecting magnification bias and the effects of 
baryons. We see immediately several features: the ef- 
fect of magnification is largest at the low £ values (small 
£ correspond to large angular scales), the magnitude of 
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the magnification correction increases with redshift, and 
the effect at high t near the baryon wiggles is mostly 
(though not completely) to boost the amplitude rather 
than change the shape of the spectrum. 

Magnification clearly changes the shape of the angu- 
lar power spectrum. The most significant changes are 
at I < 100 where the peak in the power spectrum re- 
sides. The peak in the power spectrum is an important 
feature because its location is related to the Hubble scale 
when the universe transitions from radiation domination 
to matter domination. This feature in the power spec- 
trum can be used as a standard ruler |2|. 

We identify the location £ pea k and amplitude of the 
peak in C 99 and C™". The shift in £ pea k 7 and the frac- 
tional change in amplitude, (C? n -Cf 9 )/Cf 9 (where C f n ™ 
is evaluated at the l pea k for C™ n and C 99 is evaluated at 
the £ P eak found for C 99 ) when magnification is included 
are shown in Table HI The change in £ pea k and the change 
in amplitude increase both with redshift and with in- 
creasing width of the selection function a. The value of 
(-peak is also dependent on the width of the selection func- 
tion, with £ pe ak taking slightly larger values for narrower 
selection functions. 

One additional consequence of magnification bias is to 
change the redshift and scale at which non-linear evo- 
lution of the power spectrum becomes important. This 
is because even for sources at high redshift where non- 
linearity is small, will depend on structure at low 
redshift where non-linearity is important. Indeed, even at 
zq = 3.5 the difference between Cf n calculated with the 
non-linear power spectrum and C? n calculated with the 
linear power spectrum is > 10% by £ ~ 1800, 1400, 1000 
for a = 0.07, 0.15 and 0.30 respectively. This is to be 
compared with C 99 , the no magnification case, for which 
the difference between C 99 (zq = 3.5) calculated with and 
without non-linear evolution does not approach 10% un- 
til I ~ 2600 for each value of a. 



V. THE ANGULAR CORRELATION 
FUNCTION AND THE BARYON BUMP 

We now turn our attention to the real space angular 
correlation function w nn (Q, Zq) (Equation ITT]) . In mul- 
tipole space, the baryon oscillations appear as a series 
of peaks in the power spectrum (e.g. Figure [5]). In real 
space the signature of baryon oscillations is a single bump 
in the correlation function. The location of the bump in 
w nn (9, Zq) is determined by the comoving sound horizon 
at recombination and the distance to zq. The sound hori- 
zon at recombination is measured quite precisely from the 
cosmic microwave background anisotropy, thus a mea- 
surement of 9b ao at zq can give a measure of the co- 
moving distance to zq. 

Before we address the lensing corrections to the baryon 
bump, we will first discuss a few issues that arise 
(whether or not lensing is included) when using the an- 
gular correlation function to measure the baryon oscilla- 



tion scale. The angular correlation function is in some 
sense averaging the matter power spectrum across differ- 
ent redshifts. Thus, if too broad of a selection function is 
used the baryon feature will be washed out. For a fixed 
width in redshift, the washing out is more severe at low 
redshifts because Ax ~ Az/H(z) and H(z) decreases 
with decreasing z. Additionally, the baryon bump is a 
subtle feature in the angular correlation function. It is 
customary to multiply the angular correlation function 
by 6 2 to make the baryon bump easier to identify and 
characterize. Indeed for the cases we have considered this 
is necessary in order for there to be a local maximum at 
the acoustic scale at all. Figure [3] illustrates this point: 
while there is a baryon feature, there is no local maxi- 
mum in w(9). In Figure [JJ we show 9 w(9) in the region 
near the baryon bump, in this case the baryon feature is 
quite visible for a = 0.07 but is very hard to identify for 
a = 0.30, fortunately the baryon bump becomes more 
prominent with increasing redshift. 

Since the baryon bump is broadened at low redshifts 
we limit our analysis to zq > 1.0 for a = 0.07 and 0.15 
and zq > 2.0 for a = 0.30. We have experimented with 
multiplying by 9 3 : in this case even for the broadest se- 
lection function a — 0.30 the baryon bump is visible at all 
redshifts we consider. Interestingly, the magnitude of the 
shift in the baryon oscillation scale due to magnification 
bias is sensitive to whether one looks at 9 2 w(9) or 9 3 w(9). 
The shift is much larger in the latter case, suggesting that 
the importance of magnification bias for baryon oscilla- 
tion measurements depends on precisely how the baryon 
bump is fitted. In this paper, we focus on the effects on 
9 2 w(9). 

We calculate the two-point function w nn (9, zq) by per- 
forming the sum in Equation 1121 In Figure [3] we plot 
Wnn(9) and w gg (9) with redshift bins centered at zq = 1.5 
and zq = 3.0 for a small angular range near the baryon 
bump. One can see from Figure [3] that, as expected, the 
angular correlation function is changed by magnification 
bias. 

The lensing correction to the angular correlation func- 
tion is scale dependent. To address how the baryon bump 
is changed by lensing magnification we consider the lo- 
cation of the baryon peak 9b ao j the height at the peak 
9 2 w(9bao), and the peak width. The peak width is de- 



fined to be 



(^9 2 W ) /(9 2 w) 



-1/2 



evaluated at 9bao- 



Figure [5] summarizes the fractional changes in the peak 
location, height and width as a function of (5s — 2)/b for 
selection functions with a = 0.07, 0.15 and 0.30, and 
for several redshifts. The magnitude of the changes to 
the peak increase with redshift for all a. The changes 
in peak location and width are largest at the highest 
redshifts (zq = 3 — 3.5) but the exact redshift dependence 
varies depending on the selection function. This is in 
part because broader selection functions grant a larger 
w gfi (0) term which dominates at low redshifts. In the 
region near the baryon bump, this term is more strongly 
scale dependent than w ^{6) . Additionally, the baryon 
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TABLE I: A table of quantities related to the matter-radiation equality peak for a variety of mean redshifts Zo and selection 
function widths a. The peak location (of Cf 9 ) is £ pea k, the shift in the peak location when magnification is included is 

A£ peah = l™ ak -l% ak , and the change in height at the peak is AC lpca JC lp ^ k = (c^ ak - C"f p 9 ofc ) /Cf 9 ^. For z = 0.5, we 

do not consider the case of a wide selection function a = 0.30. 
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FIG. 3: Top and middle panels: the two-point function di- 
vided by the galaxy bias w(8)/b 2 . This is using a selection 
function of width a — 0.15 centered at zo = 1.5 (top) and o = 
3.0 (middle). The dashed line neglects magnification, the solid 
line includes magnification with (5s — 2)/b = 1, the dot-dashed 
line includes magnification with (5s — 2)/b = 2. Bottom 
panel: the magnification correction, (w nn (9) — w gg (8)) /b 2 . 
The solid line is for (5s — 2)/6 = 1, the dot-dashed line is for 
(5s-2)/6 = 2 



peak is sharper and the amplitude of w gg {9) is larger for 
narrower selection functions (e.g. Figure[|]), both of these 
factors make it more difficult for magnification to effect 
the peak location. 

Magnification shifts Bbao by < 0.5% for all a and Zq 
values we considered. The changes to the width and am- 
plitude are much larger. We should point out that as the 



FIG. 4: Here we zoom in on baryon bump in 2 w(0) /b 2 at zo = 
1.5. The dashed lines neglect magnification 6 2 w gg (6) /b 2 the 
solid lines include magnification 8 2 w nn (6) /b 2 . The broader 
the selection function, the harder it is to locate the baryon 
peak, the presence of magnification bias doesn't change this. 



width of the selection function is increased the baryon 
bump is also broadened making the peak harder to iden- 
tify - this is true regardless of whether magnification bias 
is present. This is illustrated in Figure [H where we show 
9 2 w(6) for a narrow angular range for each a. One can 
see that for a = 0.30 there is no maximum at this red- 
shift. 

We have shown that the change in location of the 
baryon bump found in our calculated w nn (9) and w gg (9) 
is likely to be small. However, an actual measurement 
of 9b AO will involve fitting a curve to data points with 
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FIG. 5: Effects of magnification bias on the baryon bump 
as found in 9 2 w(8). Shown here for selection functions with 
a = 0.07 (black lowest curves), a — 0.15 (cyan middle curves), 
and a — 0.30 (red upper curves) as a function of the galaxy 
population-dependent ratio (5s — 2)/fe for several different 
mean redshifts. Upper panel: the fractional change in the 
angular location of the baryon bump (0 B n AO ~ ® 9 bao) I^bao- 



Middle panel: the fractional change in the height of the 
baryon bump, (w nn (6 B n AO ) - w gg {8 9 B 9 AO ))/w gg {6 9 B 9 AO ). Bot- 
tom panel: the fractional change of the width of the baryon 
bump - see text for how the width is defined. 



error bars. From this perspective, magnification bias - 
which adds a scale and source population dependent cor- 
rection to w gg (9) - may be of more concern. This fact is 
illustrated in the lower panel of Figure [3] Here we show 
the difference between the two-point functions with and 



without magnification (u>„ 



j). Although this dif- 



ference is slowly varying in the region around the baryon 
peak, the value of the difference is not constant. In the 
regime very near the peak the offset between w nn (9) and 
w gg {9) can be approximated as a straight line. The pre- 
cise slope of this line will depend on the population of 
galaxies via b and s. Interestingly, the magnification cor- 
rection to 9 2 w(9) is closer to a constant than it is for 
w(9). When the baryon bump in 9 3 w(9) is considered we 
find the magnification correction is more strongly scale 
dependent than for 9 2 w(9). Thus the changes to the 
baryon peak location in 9 3 w(9) are significantly larger 
(by a factor of 5 — 6). We emphasize then, that while 
Figure Ogives an indication of the effect of magnification 
bias on the baryon bump, the actual effect of magnifica- 
tion on measurements of the acoustic scale will depend 
on the method by which the peak location is fitted from 
the data. When using the BAO peak as a standard ruler 



this scale and population dependent bias should be taken 
into account. 



VI. DISCUSSION 

The angular two-point correlation function and the an- 
gular power spectrum are important cosmological statis- 
tics, but a complete understanding of systematics is nec- 
essary to use these to derive precise constraints on cosmo- 
logical parameters. Magnification bias, a gravitational 
lensing correction to the observed number density of 
sources, has been known to change the correlation func- 
tion at high redshifts [HI, [H [H 0] . We have examined 
and quantified how magnification bias changes the shape 
of the angular correlation function and power spectrum. 
The effect of magnification is to adjust both the scale- 
dependence and amplitude of the correlation function. 
We have shown that the scale-dependence changes the 
location of the matter-radiation equality peak in the an- 
gular power spectrum, and can potentially lead to a bias 
in determining the location of the baryon bump in the 
angular correlation function (see Figures [2J [31 [5] and Ta- 
ble [J). Magnification bias becomes important at high 
redshifts (z > 1.5) with the most drastic scale dependent 
changes occurring at low multipoles i < 100. Precisely 
how, if ignored, magnification would bias measurements 
of the matter-radiation equality scale or the BAO size 
will likely depend on how these quantities are extracted 
from the correlation function as well as the population of 
galaxies or quasars used for the measurement. Finally, 
the magnification terms are sensitive to non-linear evo- 
lution of structure at low redshifts even for sources at 
high redshift. Consequently, the presence of magnifica- 
tion bias changes the redshift and scale for which non- 
linear evolution of the power spectrum is visible. 

Magnification bias is a source-population dependent ef- 
fect, that is the magnitude (and even the sign) of the cor- 
rection depend on the population of galaxies or quasars 
being observed. On the one hand, one may be able to 
select a population of sources with number count slope 
s = 2/5, in which case magnification bias vanishes. How- 
ever, such a selection may reduce the number of sources 
available for analysis so this is not necessarily the best 
option. Current measurements of the galaxy angular cor- 
relation are at redshifts z < 1, so the correction from 
magnification to these observations is expected to be neg- 
ligible. Projections for measurements of galaxy cluster- 
ing from future high redshift galaxy surveys will need to 
address the effect of magnification bias. 

Finally, let us reiterate that this paper focuses exclu- 
sively on the angular correlation. The effect of magni- 
fication bias on the 3D correlation has some surprising 
new features. For instance, in certain situations, mag- 
nification bias can be important even for low redshift 
measurements. This was originally noted b y 1301 1 and is 
further analyzed in two separate papers (3ll . [32 
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